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We discuss Fermi observations of 7-ray loud binaries. We show that within hadronic model of 
activity of LSI +61° 303, detection of cut-off in the GeV spectrum constrains maximal energy of 
the primary protons. In this way, the GeV 7-ray data impose constraint on the expected neutrino 
signal (spectrum and lightcurve) from the source. We also briefly discuss perspectives of GeV band 
detection of PSR B1259-63 during the 2010 periastron passage. 



INTRODUCTION 



II. GEV SPECTRAL CUT-OFF IN 
LSI +61° 303 
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Gamma-ray-loud binary systems (GRLB) are X- 
ray binaries which emit very-high energy (VHE) 7- 
rays. Four such systems PSR B1259-63, LS 5039, 
LSI +61° 303 and HESS J0632+057, have been firmly 
detected as persistent or regularly variable TeV 7-ray 
emitters [H, 13, S 0, IE] ■ Most of the variable and tran- 
sient Galactic sources of GeV 7-rays are expected to 
belong to the GRLB class. 

The source of the high-energy activity of GRLBs 
is uncertain. It can be either accretion onto or 
dissipation of rotation energy of the compact ob- 
ject. It is commonly assumed that the 7-ray emis- 
sion is produced in result of interaction of the rel- 
ativistic outflow from the compact object with the 
non-relativistic wind and radiation field of compan- 
ion massive star. Neither the nature of the compact 
object (black hole or neutron star?) nor the geom- 
etry (isotropic or anisotropic?) and physical prop- 
erties (electron-positron or electron-proton composi- 
tion?) of relativistic wind from this compact object 
are known in in most of the GRLBs. The only excep- 
tion is PSR B1259-63 system in which the compact 
object is known to be a young rotation powered pulsar 
which produces relativistic pulsar wind. 

Two GRLBs, LSI +61° 303 and LS 5039 were found 
to be strong sources of 0.1-10 GeV 7-rays In this 

contribution we discuss constraints on the physical 
model of one of these sources, LSI +61° 303, imposed 
by Fermi observations. We demonstrate that Fermi 
observations significantly constrain the expected neu- 
trino signal from the source, within hadronic model 
of source activity. For PSR B1259-63, which is not 
(yet) detected by Fermi, we discuss the expected GeV 
emission from the source during the next passage of 
the young pulsar close to the massive star, due in 2010. 



LAT spectrum of LSI +61° 303 in sub-GeV do- 
main is described by a power law dN^/dE ~ 
with a photon index F = —2.4 6]. This slope is con- 
sistent with low-energy extrapolation of the MAGIC 
and VERITAS spectra see Fig. [TJ LAT reveals 

a suppression or cut-off in the spectrum above ^ 6 
GeV [6|. In principle, suppression of the source flux 
in 10 GeV - 1 TeV band is expected, because 7-ray 
flux in this energy band is strongly affected by the ef- 
fect of pair production on the ultra-violet photon field 
of companion star ^J- Fig[l] shows the attenuation of 
LSI +61° 303 spectrum due to the 77 pair production 
at different orbital phases under the assumption that 
the TeV flux comes directly from the compact object. 
Strongest attenuation of the source flux occurs close 
to the periastron / inferior conjunction of the orbit. 

It is clear from Fig. [T]that if the observed GeV- TeV 
emission originates from the vicinity of the compact 
object, the pair production optical depth along the 
line of sight is too small to explain the suppression 
of the source flux at several GeV. If 7-ray emission 
comes from an extended region or interaction of rel- 
ativistic wind /jet with the wind of companion star, 
the emission region could be situated closer to the Be 
star surface leading to stronger attenuation of 7-ray 
flux. Although account of spatial geometry of 7-ray 
emission region could increase the suppression of the 
source flux above ^ 10 GeV, it could hardly affect 
the flux at the energies below 10 GeV. Indeed, 7-rays 
with energies < 10 GeV could interact only with 
soft photons with energy euv ^ 25 [Ej/10 GeV] eV, 
which is much larger than the typical energy e* ~ 
3fcr* ~ 7 [T*/3 X 10* K] eV of the thermal steUar 
photons with temperature T*. 

Simple "powerlaw modified by absorption" model 
of the source spectrum is in contradiction with the 
observed orbital modulation of the signal in the 0.1- 
1 TeV energy band. Fermi orbital lightcurve exhibits 
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FIG. 1: Attenuation of LSI +61° 303 spectrum due to the 
77 pair production at different orbital phases. Blue line 
is an extrapolation of a GeV spectrum before the cut-off. 
Black, and red lines shows the effect of 7-ray absorption 
at inferior and superior conjunctions, while green and ma- 
genta lines correspond to intermediate (<^= 0.6 and 0.8) 
phases. 
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a pronounced maximum around (j) — 0-4- Fig. [T] 
shows that the effect of 77 pair production at this or- 
bital phase is small, so that a maximum at the phase 
(j) ~ 0.3 — 0.4 is expected also in the 0.1-1 TeV band 
lightcurve. To the contrary, MAGIC and VERITAS 
observe the maximum of emission at the orbital phase 
0.6 < (j) < 0-8, at which the signal should be sup- 
pressed 

Since the effect of 77 pair production can not ex- 
plain the cut-off at GeV energies and the difference 
in the orbital modulation of the source flux at GeV 
and TeV energies, one has to consider a possibility 
that GeV and TeV 7-rays are produced via different 
mechanisms (e.g. synchrotron and inverse Gompton) 
and/or by different particle populations (e.g. elec- 
trons and protons). Fig. [51 adopted from the Ref. 
p^ . shows two possible models of broad-band spec- 
trum, consistent with Fermi data. In the upper panel 
the GeV band emission is inverse Gompton emission 
from electron distribution with a high-energy cut-off 
at Ecut — 10 GeV. Alternatively, in the model shown 
in the lower panel of Fig. [21 the GeV band emission 
is synchrotron emission from electron distribution ex- 
tending to the energies ^ 100 GeV (see [Toj). 

In the synchrotron model of GeV emission, shown in 
the lower panel of Fig. [2l the electron injection spec- 
trum has to be very hard 0, [l3| . Such hard electron 
spectrum could be produced if electrons are directly 
injected from relativistic pulsar wind or produced 
in interactions of high-energy protons jSj . 



III. ORBITAL MODULATION OF 
NEUTRINO FLUX 

If high-energy electrons emitting GeV synchrotron 
radiation are injected via interactions of high- 
energy protons, measurement of a cut-off energy 



FIG. 2: Model spectra compared to the data, which are the 
same as those in [3| except for the added data from the 
VERITAS ^ and Fermi Q] telescopes. The dot-dashed, 
dotted and dashed curves show the spectral components 
from the synchrotron, bremsstrahlung and IC processes, 
respectively. Adopted from [lol |. 
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FIG. 3: Illustrative scheme for the restricted time of neu- 
trino emission. Grey shaded region shows the "shadow" of 
the Be star disk on the orbital plane, from the viewpoint 
of observer. 



of Fermi spectrum E-y^cut — 6 GeV 6] imposes 
a constraint on the high-energy cut-off in the 
spectrum of the secondary e^e^ pairs, -Be, cut — 
10^^ [B/W G]~'^^'^ [E-y^aut/G GeY]^^^ eV, where B is 
the magnetic field strength in the emission region. 
This, in turn, constrains the cut-off in the spec- 
trum of primary protons, Ep^cut — K~^Ee^cut — 
0.5[B/10 G]"^/^[£:.^,cut/6 GeV]^/^ PeV where k ~ 
0.2 is the typical inelasticity of pp collisions. 

PeV protons do not interact efficiently with the UV 
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photons from the Be star, since the center of mass en- 
ergy of proton-photon colhsions is below the thresh- 
old of pion production Ep^ ~ 2 x 10^^ [e*/10 eV] eV. 
Therefore Fermi data favor proton-proton hadronic 
models as compared to proton-photon interaction 
models, unless the magnetic field in the interaction 
region is at the level of i? < 0.1 G. 

The synchrotron interpretation of the source spec- 
trum in the keV-GeV energy range enables to de- 
duce the shape of the electron spectrum from the 
measured shape of the synchrotron spectrum. This 
was done in the Ref. Q in which it was found 
that the data are consistent with the assumption 
that e~^e~ pairs are initially injected with energies 
close to the high-energy cut-off. Hard injection spec- 
trum of e+e" pairs could be explained by the ini- 
tial hard (much harder than E^'^) injection spec- 
trum of protons (e.g. injection from the proton- 
loaded relativistic wind with large bulk Lorentz fac- 
tor). The spectrum of neutrinos produced in in- 
teractions of high-energy protons is expected to be 
sharply peaked at the highest energies E^, ~ KEp ~ 

IQi'* [B/10 G]"^/^ [E^,cut/& GeV]^/^ TeV. 

Neutrino flux is expected to be modulated on the 
orbital time scale. Accounting for the orbital period- 
icity of the neutrino signal significantly increases the 
IceCube potential for the detection of the source, once 
the validity of the model prediction is assumed. The 
resulting model-dependent approach in search for a 
neutrino signal is diametrically opposed to a random 
search for steady neutrino point source over the whole 
sky and complementary to a search for steady neutrino 
point source selected out of a catalogue [11| . 

The neutrino flux is not expected to be positively 
correlated with the TeV emission, which is affected 
by the orbital phase dependent absorption in inter- 
actions with UV photons from Be star (see above). 
It is also not expected to be correlated with the GcV 
band 7-ray emission, which might be strongly affected 
by development of electromagnetic cascade initiated 
by high-energy proton interactions. Thus, neutrino 
lightcurve has to be calculated independently, rather 
than read out from photon lightcurve. 

PeV protons are able to penetrate deep into 
(or even through) the dense equatorial disk of Be 
star, because their Larmor radius, Rl ~ 3 x 
10^2 [Ep/10^^ eV] [B/1 G]~^ cm, is comparable to 
the thickness of the disk. Since proton trajectories 
are not randomized by the magnetic fields in the disk, 
the neutrino emission from the system is expected to 
be anisotropic, with the maximum flux emitted into 
a cone with the axis aligned with the compact object 
- Be star direction [2]. A neutrino signal could be 
detected if the neutrino emission cone passes through 
the line of sight, i.e. if the compact object is situated 
behind the dense equatorial disk of Be star during a 
part of the orbit, as it is shown in Fig. [3l 

We consider a disk which is truncated at the dis- 
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FIG. 4: The thick lines represent the orbital phase of the 
superior conjunction and periastron. The thin curves show 
the range of orbital phase with an expected neutrino emis- 
sion for different sizes of equatorial disk around Be star. 



tance i?disk comparable to the binary separation of the 
star and compact object at the periastron The 
span of the orbital phase period during which the com- 
pact object is situated behind the disk is maximized 
if the disk is oriented perpendicularly to the line of 
sight, as it is shown in Fig. [3] For such a disk orienta- 
tion, the result of calculation of the interval of orbital 
phases, during which the compact object is in the disk 
shadow (i.e., the interval during which neutrino signal 
from the source could be detected), is shown in Fig. 
[31 For this calculation, we use the estimates of the 
orbital parameters of the system from Ref. [l^l ■ 

The interval of phases at which neutrino signal is ex- 
pected, depends on the (unknown) inclination of the 
binary orbit. The width of the phase interval of neu- 
trino emission increases with the increase of the disk 
size. In Fig. 0] we show the results of calculation of 
the phase intervals of neutrino emission for a range 
of disk sizes, from i?disk equal to the half of the bi- 
nary separation distance at periastron Dpcr to twice 
the binary separation at periastron. 

It is interesting to note, that the equatorial disk 
of Be star in the LSI -f61° 303 system is observed 
to periodically vary on super-orbital time scale of 
Tso — 4.6 yr [1J|. This means that the orbital phase 
interval of neutrino emission is also expected to be 
modulated on the 4.6 yr time scale. The strongest 
neutrino signal (corresponding to the largest span 
of the neutrino emission phase interval) is expected 
when the equatorial disk expands to its maximal size, 
roughly once in each 4.6 yr cycle. 

It is clear from Fig. 0] that, independently of the 
disk size and inclination of the orbit, the phase inter- 
val of neutrino emission from the source is 6(j) ^ 0.3. 
If this neutrino emission model within a reduced time 
window closely matches the actual source neutrino 
emission, the sensitivity of a neutrino telescope for 
the detection of the source is increased due to (1) an 
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atmospheric neutrino reduction by a factor 54>^^ ^ 3 
and (2) beyond a free phase variabihty fit due to the 
constraint of the phase of the neutrino emission. 

IV. ORIGIN OF THE PSR B1259-63 
EMISSION. 

PSR B 1259-63 is the only 7-ray loud binary sys- 
tem for which the nature of the compact object and 
the source of the power are certain. It is clear that 
(unpulsed) radio-to- TeV emission from the system is 
generated in collision of relativistic pulsar wind with 
the dense equatorial wind of Be star. The line of inter- 
section of the disk plane and the orbital plane is tilted 
at about 90° with respect to the major axis of the bi- 
nary orbit, so that the pulsar passes through the disk 
twice per orbit. Episodes of pulsar passage through 
the disk are associated with pronounced flares of the 
source, visible in radio, X-ray and TeV 7-ray bands 

[i|,[il[iE[il,E3- 

Similarly to the case of LSI -1-61 303, physical mech- 
anism through which X-ray to 7-ray emission from the 
source is produced (synchrotron or inverse Compton), 
is not well constrained [3, [H, [IB, [3 ■ In fact, spec- 
tral and timing characteristics of the source could be 
successfully reproduced in both models, see Fig. [51 

From Fig. [5] it is clear that Fermi observations of 
the next flaring period of the source during 2010 pe- 
riastron passage will clarify the mechanism of X-ray- 
to- 7-ray emission. Predictions of the two models sig- 
nificantly differ in the two cases. If the X-ray-7-ray 
emission is dominated by inverse Compton emission 
from single population of electrons, the source should 
be detected by Fermi during the several month-long 
flare around the periastron. The source spectrum is 
expected to be a powerlaw. At the same time, the 
source could remain below Fermi sensitivity level if 
the the X-ray emission is dominated by synchrotron. 
The cut-off energy of the synchrotron spectrum is not 
currently constrained by the available data. If the en- 
ergies of electrons accelerated at the interface of pulsar 
and stellar wind (or injected directly from the pulsar 
wind) reach ~ 100 TeV energies (similarly to the case 
of LSI -1-61° 303), the high-energy cut-off of the syn- 



chrotron spectrum might be revealed by Fermi. 



V. SUMMARY 

We discussed the results of Fermi observations of 
GRLBs. We found that Fermi data significantly 
constrain the mechanisms of 7-ray emission from 
LSI -f61° 303 and enable to work out firm predic- 
tions for the expected neutrino signal from the source 
in the framework of hadronic models of source activ- 
ity. For PSR B1259-63 we discussed the possibility 
of detection of the GcV flare from the source in 2010, 
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FIG. 5: Spectral energy distribution of PSR B1259-63 in 
the IC (top) and synchrotron (bottom) models (from Ref. 
[l3l). Thick cyan line shows estimated 1 month Fermi 
sensitivity. 



during the period of passage of the young pulsar in 
the system close to the massive star. 
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